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ABSTRACT 


,,  The  effect  of  finite  chemical  reaction  rates  on  the  chemical  structure 
of  hypersonic  wakes  is  examined.  The  technique  of  incipient  non-equilibrium 
coefficients,  based  on  linearized  chemical  kinetics  for  a 


component 


(Og,  Ng,  0,  N,  NO,  e;;j)  tw  reaction  model  of  high  temperature  air,  is 


utilized  in  the  analysis.— 

^.C^The  far  wake  is  treated  as  an  asymptotic  region;  the  near  wake  theory 
considers  the  viscous  mixing  of  the  boundary  layer  flow  with  the  inviscid  and 
re-circulating  flows  in  the  base  region,  and  the  subsequent  development  of 
this  viscous  core  into  the  wake  region  downstream.  A  range  of  altitudes  and 
flight  speeds  of  practical  interest  is  considered.  <  The  results  show  where 
non- equilibrium  effects  are  important  in  the  wake  and  whether  they  arise 
primarily  from  convective  or  diffusive  effects.  The  species  that  are  the 


principal  contributors' to  the  non-equilibrium  effects  are  also  identified.-- 


'■5 


The  incipient  non-equilibrium  results  are  then  related  to  calculations 
assuming  chemically  frozen  flow. (/  By  the  application  of  chemical  kinetics 
to  the  frozen  results,  it  can  be  shown  that  the  appearance  of  incipient  non¬ 
equilibrium  effects  bears  an  inverse  relationship  to  the  validity  of  the  frozen 
flow  assumption  for  equivalent  wake  calculations. 


NOTE:  This  work  was  performed  under  Contract  Number  AF  04(647) -6 17. 
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INTRODUCTION 


Aerodynamics  at  hypersonic  speeds  is  greatly  complicated  by  chemical 
changes  occurring  in  high  temperature  air.  These  chemical  changes  have 
significant  effects  on  the  thermodynamic  and  transport  properties  of  air; 
in  some  cases,  they  give  rise  to  diffusive  phenomena  which  can  contribute 
to  the  transfer  of  energy  through  the  moving  and  chemically  reacting  fluid. 

In  addition,  the  theoretical  prediction  of  the  chemical  state  of  air  in  the 
wake  is  an  important  task  for  its  own  sake,  because  it  determines  the 
electrical  transmission  and  radiative  properties  of  the  wake,  which  are 
important  means  for  the  detection,  observation,  and  identification  of 
hypersonic  re-entry  bodies. 

For  these  reasons  a  great  deal  of  current  theoretical  aerodynamic 
research  has  been  directed  towards  predicting  the  chemical  state  of 
reacting  air  flows.  Under  many  circumstances  of  practical  interest  in 
the  re-entry  problem,  chemical  changes  will  take  place  at  rates  which  are 
comparable  to  rates  at  which  changes  occur  in  the  flow  due  to  fluid  mechanic 
phenomena.  Under  these  circumstances  theoretical  predictions  must  be 
based  on  simultaneous  consideration  of  chemical  kinetics  and  fluid  mechanics. 
This  is  a  very  formidable  problem,  particularly  in  diffusive  flows  and  for  a 
chemical  kinetic  model  that  attempts  to  account  for  all  the  ionizing  species  in 
a  physically  meaningful  way.  For  these  reasons  much  of  current  research 
has  been  directed  towards  justifying  simplifications  that  make  the  general 
problem  of  coupled  chemical  kinetics  and  fluid  mechanics  more  tractable  in 
specific  situations. 

The  present  investigation  is  directed  toward  solving  this  general  problem 
for  the  specific  case  of  uncontaminated  air  wakes  behind  bodies  moving  at 
hypersonic  speeds.  This  type  of  flow  is  characterized  by  the  presence  of 
both  convective  and  diffusive  effects,  but  the  latter  (as  in  boundary  layers) 
occur  only  along  constant  pressure  surfaces.  An  approximation  procedure, 
which  is  based  on  expanding  the  coupled  fluid  mechanic  and  chemical  kinetic 
equations  in  powers  of  a  small  aerothermochemical  parameter,  is  used  to 


describe  these  types  of  flow.  The  first  term  in  the  expansion  describes 
the  hypersonic  wake  corresponding  to  instantaneous  local  chemical 
equilibrium.  The  next  term  in  the  expansion,  which  is  calculated  in 
detail  for  several  practically  significant  wake  flows,  indicates  the  deviation 
of  the  flow  from  the  chemical  equilibrium  model  due  to  the  finite  rate 
chemical  kinetics.  These  deviations  are  due  both  to  the  rapidity  of 
convection  of  the  fluid  particles  along  the  wake  streamlines  and  the  inter¬ 
diffusion  of  the  reacting  species  across  the  streamlines.  The  results  show 
which  of  these  non -equilibrium  effects  dominates  in  various  parts  of  the 
wake  and  also  which  of  the  species  are  the  major  source  of  these  effects. 
The  analysis  is  based  on  a  seven  component  (Og,  Ng,  O,  N,  NO,  NO'*’,  e') 
ten  reaction  model  of  high  temperature  air. 

METHOD  OF  ANALYSIS 

The  key  equation  that  couples  the  chemical  kinetics  and  the  fluid 
mechanics  in  reacting  flow  problems  is  the  equation  of  conservation  for 
the  chemical  species  that  comprise  the  reacting  mixture.  For  a  mixture 
comprised  of  n  +  m  species  with  m  independent  conservation  laws  relating 
the  species  concentrations  to  each  other  (e.  g.  over-all  mass  conservation, 
charge  conservation  in  a  neutral  plasma,  etc. )  there  will  be  n  independent 
species  conservation  laws  of  the  following  type; 


DC .  .  w, 

— '  +  -divpC.i.  =  -^(p.T.C,  ••) 


(1) 


In  addition  to  equations  (1)  the  usual  fluid  mechanic  conservation  laws 
for  mass,  momentum,  and  energy,  the  equation  of  state  and  other  auxiliary 
equations  relating  thermodynamic  a;id  transport  properties,  are  necessary 
to  complete  formulation  of  the  general  problem.  However,  due  to  the 
appearance  of  the  chemical  kinetic  mass  generation  rate  terms  on  the  right 
hand  side  of  equations  (1)  these  equations  are  of  particular  significance  in 
the  following  analysis. 
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In  reference  1  it  was  pointed  out  that  for  "fast"  chemical  reactions  it 

is  possible  to  treat  the  general  problem  in  terms  of  an  expansion  procedure. 

The  expansion  parameter  for  the  procedure  is  the  ratio  of  an  appropriate 

chemical  time,  t  (which  is  based  on  the  sensitivity  of  the  chemical  mass 
c 

generation  rate  functions  to  changes  in  composition)  and  a  fluid  mechanic 
reference  time  t^.  As  a  result  of  the  expansion  procedure  incipient  non¬ 
equilibrium  coefficients  can  be  defined  for  a  reacting  mixture,  which  can 
be  used  to  predict  non- equilibrium  effects  both  qualitatively  and  quantitatively 
from  properties  of  the  flow  solution  based  on  instantaneous  local  chemical 
equilibrium.  The  latter  type  solutions  are  comparatively  simple  to  obtain 
for  wakes  and  other  hypersonic  aerodynamic  problems  of  current  interest. 

In  what  follows  the  results  of  reference  1  will  be  utilized  in  conjunction 
with  equilibrium  wake  solutions;  for  proper  interpretation  of  the  results  it 
is  appropriate  now  to  reiterate  the  salient  features  of  the  expansion 
procedure  including  the  concept  of  incipient  non- equilibrium  coefficients 
as  they  apply  to  wake  type  flows. 

If  equations  (1)  are  written  in  non-dimensional  form  the  expansion 
parameter  ^ 

e  =  -p  (2) 

'f 

will  appear  explicitly  in  the  equations  as  follows; 


e 
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div(pc.5,) 


=  T 


(3) 


If  now  all  dependent  variables  are  expanded  in  powers  of  6  and  terms 
grouped  according  to  powers  of  £  in  all  the  pertinent  equations,  the  over-all 
conservation  equations  for  the  leading  terms  are  unchanged  (since  the 
parameter  £  does  not  appear  in  these  equations  explicitly),  but  species 
conservation  equations  (3)  become  for  the  leading  term 


w. 

(p,,T,,C.^)  =  0  (4) 


3 


where  subscript  1  denotes  the  leading  terms  in  the  £  expansion.  Equations 
(4)  form  an  algebraic  system  which  can  be  solved  for  the  first  order 
composition 


C. 


(Pi.  Tj ) 


(5) 


without  reference  to  any  fluid  mechanic  problem,  thereby  decoupling  the 
chemical  problem  from  the  fluid  mechanics.  This  justifies  use  of  the 
chemical  equilibrium  composition  for  these  types  of  problems.  The  effect 
of  the  flow  process  on  composition  can  be  ascertained  by  considering  the 
next  term  in  the  expansion  of  equations  (3).  Specializing  now  to  axi- 
symmetric  steady  wakes  the  next  set  of  equations  in  the  expansion  of  (3) 
become: 
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This  is  a  system  of  n  linear  algebraic  equations  for  the  unknown  second 
order  terms;  formulation  of  the  second  order  problem  should  be  completed  by 
considering  corresponding  expansions  of  the  remaining  conservation  equations. 
However,  using  (5)  in  the  left  hand  side,  and  (4)  in  the  right  hand  side  of  (6) 
these  equations  can  be  written  in  a  more  immediately  useful  form: 
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where  the  quantities 

AC.  H  C. 
J2  J2 
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group  all  second  order  terms  into  one  "composition  deviation"  term,  and 
the  terms 


Pc  ^ 


•  •  aT,  ,  aT, 

T  s  u,'  — ^  +  V, 

CD  ^  ax  1  by 


1  ^  tfiY  5T, 

ypi  ay  Sc  ay 


t  =  -  -1^ 

D  Sc  \  ay  / 


account  for  the  flow  effects  based  on  the  first  order  (i.  e.  chemical  equilibrium) 

solution  of  the  wake.  The  linear  algebraic  system  (7)  is  solvable  for  the  AC.  , 

(2) 

furthermore,  as  pointed  out  in  an  earlier  note  ,  the  quantities 


can  be  calculated  once  and  for  all  without  considering  any  particular  flow 
problem.  Associating  each  of  the  species  with  the  appropriate  amount  of 
chemical  energy,  which  is  required  to  create  it  from  a  low  temperature 
reference  state,  incipient  non -equilibrium  coefficients  can  be  defined  which 
account  for  the  combined  effect  of  all  the  species; 


AE  _  r 
R-  '  I 


n+m  AC, 

^  -f  E. 

^CD  A  ‘ 


is.  -  T  V 
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For  a  given  equilibrium  wake  solution  the  local  value  of  enthalpy  is 
known  and  the  quantities  in  (9)  can  be  calculated,  so  that  the  ratio 


AE  1  AE  .  AE  .  .  AE  ^ 

»  ■  ^  T  rn  T-  D 


can  be  calculated  everywhere  in  the  flow.  This  ratio,  AE/h,  can  be 
considered  an  indicator  of  incipient  non-equilibrium  effects  in  the  flow. 

If  it  is  small  compared  to  1.0  the  equilibrium  solution  is  valid.  If  the 
ratio  is  of  the  order  or  larger  than  1.0  non-equilibrium  effects  are  present 
in  the  flow.  Due  to  the  linearity  of  all  the  pertinent  equations  the  causes 
of  deviation  from  equilibrium  can  be  ascertained  both  as  to  the  fluid 
mechanic  effect  (i.  e.  convection,  diffusion,  etc. ),  from  the  magnitude 
of  the  various  terms  in  (9);  and  as  to  the  chemical  species  that  are  the 
principal  cause  on  the  deviation  from  equilibrium,  from  the  magnitude  of 
the  various  terms  in  (11). 

Incipient  non-equilibrium  coefficients  (11)  have  been  calculated  for  a 

seven  component  (Og,  Ng,  0,  N,  NO,  NO^  e* )  ten  reaction  model  of  high 

(31 

temperature  air  given  by  Bortner  and  Golden'  :  The  ten  reactions  that  were 
considered  are  reiterated  below  for  the  sake  of  reference: 


o  +  Og  - 

Og  +  M3 
0  +  Ng  :: 

O  +  NO  ~ 
N  +  0  ~ 


N  +  N, 


Ng  +  Mg 


20  +  O, 


20  +  M3 

N  +  NO 
N  +  0, 
NO+  +  e- 
2N  +  N, 


2N  +  Ma 


M,  =  N,  N„,  NO,  NO+ 


Mg  =  0,  0„,  NO,  NO+ 


(10)  NO+M^o 


N  +  O  +  Mjo  Mio=  0,02,N,  N2,N0+ 


The  results  of  these  calculations,  which  are  described  in  more  detail 


in  reference  1,  were  applied  to  two  different  types  of  wake  calculations. 
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WAKE  ANALYSES 


A  simple  yet  informative  description  of  the  far  wake  can  be  obtained 
by  an  analysis  based  on  the  asymptotic  behavior  of  the  disturbed  fluid  far 
from  the  body.  The  basic  equation  for  this  analysis  is  the  momentum 
conservation  integral: 

C 

where  u^^  is  the  velocity  defect  with  respect  to  the  free  stream  velocity; 

u  =  u  -  u  (15) 

D  ® 

and  y  is  the  radial  distance  from  the  centerline.  The  velocity  and  enthalpy 
are  related  by  means  of  the  adiabatic  energy  relation: 


4ff 


J 


£.  ^ 
0  u 


ydy 


(14) 


h  +  y  =  H  (16) 

Ct  00 

The  total  enthalpy  of  the  flow  H  and  the  free  stream  pressure  are  both 

00 

assumed  to  be  constants.  A  compressibility  transformation  and  Prandtl's 
mixing  length  analysis  are  then  applied  so  that  the  description  of  the  far  wake 
is  defined  by; 


where  rj  is  the  transformed  radial  coordinate; 

(18) 

^00 
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Equations  (17)  and  (14)  determine  the  variation  of  Rwith  x: 

V3 


(19) 


The  details  of  the  analysis  are  set  forth  in  reference  4.  Due  to  the  lack  of 
pressure  gradient  the  quantity  p-  will  be  zero;  the  other  two  quantities 
in  (17)  can  be  calculated  from  the  above  analysis.  (The  turbulent  Schmidt 
number  is  assumed  to  be  one;  this  is  consistent  with  the  remaining 
assumptions. ) 

In  the  near  wake  analysis,  the  integral  forms  of  the  continuity,  momentum 
and  energy  equations  have  been  used.  They  are,  when  the  external  inviscid 
velocity  and  enthalpy  gradients  are  neglected  (reference  5). 

Continuity  Equation; 

^[p  u  6^A]=p  u6  [*•—  -  —1  (20) 

dx  ^  '^e  e  ■'  "^e  e  j^dx  u^  J 

Momentum  Equation; 


_d 

dx 
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(21) 


Energy  Equation; 
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u  6®  M0,  , 

e  h  dx 
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where 
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and  6  is  the  viscous  wake  radius. 


In  addition,  momentum  and  energy  balances  along  the  axis  gave  two 
more  equations  as: 


3u 


.i  +  liM 
2  dx  I  Re 
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£u  1 


(24) 
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dx 
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Re  Pr 
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(24) 
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Non-similar  profiles  for  pu  and  h  were  then  assumed,  i.  e. , 


^  =  F(Xj7))  +  a^{x){l-rf)  +  a3(x){l-Tj®)®+a3(x)(l-Tf)^ 

"e  e 

(25) 


~  =  G(x,T,)+b,{x)(l-r)?)  +  b,(x){l-rf)"+ba(x){l.r;2)" 
e 

Equations  (20)  to  (25)  are  the  basis  of  the  near  wake  analysis.  To 
start  the  near  wake  calculations,  initial  conditions  near  the  "neck"  must 
be  known.  In  the  present  work,  these  conditions  were  obtained  by  using 
the  analysis  for  base  flow  given  in  Reference  5.  In  that  base  flow 
analysis,  the  initial  finite  boundary  layer  at  the  comer  of  the  base  and 
its  subsequent  mixing  with  the  inviscid  and  recirculating  flows  in  the 
base  region  have  been  included  in  the  treatment.  Other  factors  affecting 
the  solution  are  the  Reynolds  number,  base  wall  enthalpy,  base  pressure 
values,  and  the  inviscid  flow  conditions  that  were  assumed.  Therefore, 
these  results  give  a  more  detailed  and  realistic  description  of  the  wake 
structure  in  the  vicinity  of/and  further  downstream  of  the  "neck"  region 
than  can  be  obtained  by  means  of  the  simplified  treatment  outlined  above 
for  the  far  wake.  For  particular  near  wake  solutions  the  quantities  in  (9) 
were  then  calculated  in  order  to  assess  incipient  non- equilibrium  effects. 

DISCUSSION  OF  METHOD 

It  is  quite  evident  that  the  expansion  procedure  that  forms  the  basis  of 
the  incipient  non- equilibrium  analysis  is  not  valid  unless  succeeding  terms 
in  the  expansion  are  small  throughout  the  flow  field.  This  implies  that  if 
deviations  from  chemical  equilibrium  occurred  in  the  flow  field  surrounding 
the  re-entry  body,  application  of  this  method  to  the  wake  could  not  be 
justified  rigorously.  However,  the  sensitivity  of  the  chemical  rate  effects 
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to  pressure  and  temperature  seems  to  indicate  that  upstream  effects 
play  a  minor  role  compared  to  local  properties  of  the  flow  in  bringing 
about  non-equilibrium  effects.  This  will  tend  to  justify  the  use  of 
incipient  non-equilibrium  coefficients  on  a  local  basis;  it  might  also 
permit  identification  of  chemically  frozen  flow  regions  where  the  values 
of  the  incipient  non-equilibrium  indicators  AE/h  are  large. 

In  order  to  ascertain  whether  this  is  indeed  so  chemically  frozen 
wake  flow  was  assumed  for  the  same  flight  conditions.  The  validity  of 
this  assumption  was  then  tested  by  using  the  frozen  composition  and  an 
appropriate  local  temperature  to  calculate  the  chemical  mass  generation 
rate  terms  for  all  the  species.  These  terms  were  then  compared  to  the 
diffusive  terms  in  the  species  conservation  equation  on  the  basis  of  the 
chemical  energy  involved  in  all  the  species,  The  calculations  utilize 
the  fact  that  in  the  absence  of  pressure  gradient  and  for  a  Schmidt  number 
of  1.0  the  relation 


Fr 


=  C.  +  —  (C.  -  C.  ) 


(26) 


ref 


ref' 


identically  satisfies  the  frozen  species  conservation  differential  equation 

provided  u  ,  C.  i  and  C.  are  all  constants,  independent  of  x  and  y. 

»  ^  ^ref 

Then  an  "energy  generation  rate",  due  to  the  reactions 


n+m  w. 


E„  5  y  E.  —  (p.T.C.  \ 

1  =  1  '  * 


(27) 


was  compared  to  an  "energy  accumulation  rate";  due  to  diffusion; 
.n+m 

E-  5  y  E.  —  —  ,iy - (28) 

D  Z+  1  py  By  9y 


i=l 


•  t 


If  the  ratio  E  /E  is  small  frozen  flow  is  justifiable.  The  outcome 

IV  D 

of  this  comparison  depends  greatly  on  the  temperature  chosen  to  calculate 


11 


the  rates  in  (27)  and  on  the  reference  condition  that  is  used  to  calculate 

the  frozen  compositions  in  (26).  For  the  far  wake«  the  most  "sudden" 

freezing  would  occur  if  the  wake  composition  froze  at  a  rear  stagnation 

point  behind  the  base.  If  the  C{  are  taken  as  the  equilibrium  compo> 

ref 

sition  corresponding  to  ambient  pressure  and  the  free  stream  stagnation 
enthalpy,  the  frozen  wake  composition  profiles  would  exhibit  the  maximum 
possible  deviation  from  the  equilibrium  assumption.  Two  temperatures 
were  then  used  to  calculate  the  reaction  rates;  the  temperature  corres¬ 
ponding  to  the  equilibrium  far  wake  solution  is  the  lowest  possible 
temperature,  the  actual  temperature  might  be  higher,  the  maximum 
being  a  temperature  corresponding  to  the  completely  frozen  composition 
described  by  equation  (26).  Neglecting  any  effect  of  composition  on  the 
equation  of  state,  the  enthalpy  profiles  in  the  far  wake  were  assumed  to 
be  unchanged  from  equations  (16)  and  (17),  so  that  a  "frozen  temperature" 
could  be  calculated  by: 


n+m 


where  the  species  enthalpies  are  given  functions  of  temperature. 

Steps  (26)  through  (29)  were  performed  for  the  near  wake  also.  The 
free  stream  condition  in  (26)  was  replaced  by  the  external  (inviscid) 
condition  used  in  the  near  wake  analysis,  and  the  appropriate  reference 
conditions  needed  in  (26)  was  taken  as  the  surface  condition  at  the  trailing 
edge  of  the  body.  (This  corresponds  to  a  catalytic  wall  assumption. )  The 
remaining  assumptions  of  unchanged  enthalpy  and  velocity  profiles  (as 
compared  to  the  equilibrium  solution)  were  made  again.  The  presence 
of  a  pressure  gradient  in  the  near  wake  (which  invalidates  result  (26) ) 
was  minimized  by  an  appropriate  choice  of  external  conditions. 
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NUMERICAL  RESULTS 

Numerical  calculations  were  performed  for  the  far  and  near  wakes 
for  the  flight  conditions  listed  in  Table  I.  For  the  equilibrium  far  wake, 
velocities,  temperatures,  enthalpies,  and  chemical  energies  along  the 
centerline  are  given  in  Figures  1  through  3,  respectively.  The  chemical 
energy  (Cchem^  energy  absorbed  by  dissociation  and  ionization; 

when  compared  to  h  it  is  an  indicator  of  the  potential  significance  of 
chemical  effects.  Figure  4  shows  ^E/h.,  the  incipient  non- equilibrium 
index  (I.  N.  I. ),  for  the  same  six  cases  abng  the  far  wake  centerline. 

Since  the  temperatures  are  the  highest  near  the  centerline,  for  a  1  foot 
body^ =  l.oj  this  figure  indicates  the  existence  of  an  equilibrium 
region  with  the  far  wake  (where  the  I.  N.  I.  is  less  than  1.0)  at  only  3  of 
the  6  conditions  studied.  Near  the  body  N  is  the  largest  contributing 
species  to  the  I.  N.  I. ,  whereas  further  downstream  it  is  0.  It  is  note¬ 
worthy  that  as  the  temperature  drops,  the  I,  N,  I.  first  takes  a  downward 
turn  for  many  of  the  cases  and  then  begins  to  rise.  This  phenomenon  is 
due  to  the  fact  that  in  spite  of  the  lower  rates  at  the  lower  temperatures, 
the  tendency  to  deviate  from  equilibrium  becomes  smaller  because  of  the 
diminished  supply  of  N  atoms  that  are  "required"  to  recombine,  and  the 
simultaneous  "flattening"  of  the  wake  profile  which  implies  a  decrease 
in  the  rates  that  are  "prescribed"  by  the  equilibrium  wake  solution  (i.  e. , 
the  terms  in  equation  (9)  become  smaller).  However,  as  soon  as  O 
recombination  begins  to  take  place,  the  tendency  to  deviate  from  equilibrium 
becomes  several  orders  of  magnitude  larger.  One  may,  for  example, 
surmise  from  these  results  that  for  the  high  altitude  condition  (case  1), 
both  the  N  and  O  recombinations  will  lag  far  behind,  whereas  at  lower 
altitudes  (e.  g. ,  cases  3  and  4)  most  of  the  N  atoms  will  have  recombined 
along  the  centerline;  the  excess  atoms  in  the  far  wake  will  be  mostly  0 
atoms.  This  might  have  important  consequences  regarding  electron 
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densities,  due  to  the  increased  occurrence  of  NO'*’  ions  in  regions  where 
both  type  of  atoms  are  present. 

Figures  5  through  7  show  properties  of  a  particular  equilibrium  far 
wake  profile  as  to  velocity,  temperature,  enthalpy,  chemical  energy  and 
composition.  A  station  383  downstream  from  the  body  for  case  3 

was  chosen  as  the  example.  Figure  9  shows  the  corresponding  I.  N.  I. 
profile,  showing  also  the  separate  effects  of  convective  and  diffusive 
terms.  It  is  clear  that,  except  at  the  wake  center,  interdiffusion  of  the 
reacting  species  is  the  dominant  reason  for  deviating  from  chemical 


equilibrium  in  the  far  wake.  This  inter  diffusion  is  due  both  to  the 

(a®  T\ 

j ,  and  to  the  curvature  of 

the  equilibrium  composition  vs.  temperature  variation  in  conjunction  with 
the  temperature  gradients  as  shown  in  Figure  10.  (Curvatures  of 

the  equilibrium  composition  profiles  (Figure  7)  also  indicate  the  regions 


where  there  will  be  no  diffusion. )  The  breakdown  of  the  convective  effect 


into  axial  and  radial  components  (Figure  11)  further  demonstrates  the 
utility  of  the  method.  Figure  12  indicates  contribution  of  the  various 
species  to  the  I.  N.  I. ,  atomic  N  and  O  dominate  near  the  center  and  in 
the  outer  region  of  the  far  wake,  respectively;  the  contribution  of  NO 
appears  to  be  significant  only  at  a  narrow  intermediate  region. 

Figure  8  compares  I.  N.  1.  profiles  in  the  far  wake  for  three  different 
cases.  Comparison  of  cases  3  and  4  shows  the  effect  of  altitude;  whereas 
comparing  cases  4  and  5  the  effect  of  flight  speed  at  the  same  altitude  can 
be  ascertained.  As  expected,  low  altitudes  and  high  speeds  (and  the 
ensuing  high  temperatures)  will  extend  the  regions  of  equilibrium  flow 
both  laterally  and  axially.  This  is  even  more  clearly  shown  on  Figure  13 
where  the  contours  of  ^E/h  =  1.0  are  plotted  on  a  wake  map  for  the  three 
cases  where  equilibrium  regions  did  occur.  The  effect  of  body  size  on  these 
contours  for  one  of  the  cases  is  also  shown,  indicating  that  the  equilibrium 
flow  regions  will  be  different  for  various  size  and  shape  re-entry  bodies. 
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Similar  calculations  were  performed  for  several  near  wake  profiles. 

Temperatures,  enthalpies  and  chemical  energies  are  shown  in  Figures  14 

and  15,  respectively.  The  I.  N.  I.  profiles  are  shown  in  Figure  16.  They 

show  that  for  the  high  altitude  condition  (cases  7  and  8)  the  equilibrium 

assumption  is  not  at  all  reasonable.  It  is  less  reasonable  than  in  the 

corresponding  far  wake  (case  2),  in  spite  of  the  considerably  higher 

temperatures  and  somewhat  higher  pressure  near  the  base.  The  reason 

for  this  is  partly  the  more  rapid  variation  of  fluid  properties  due  to  the 

large  gradients  near  the  base  (i.  e. ,  T_  and  T  are  large)  but  mainly 

CO  o 

the  fact  that  at  this  pressure  there  is  a  great  deal  of  variation  in  the  Ng  -  N 
equilibrium  at  the  higher  temperatures,  whereas  (as  mentioned  previously) 
at  the  lower  temperatures  only  small  chemical  change  is  "required"  by 
the  equilibrium  solution  until  the  temperature  is  sufficiently  low  for  the  0 
recombination  reaction  to  begin.  At  the  lower  altitude,  the  equilibrium 
assumption  appears  justified  in  some  regions  of  the  near  wake  at  high  wall 
temperatures  (case  10).  For  case  9,  Figure  17  indicates  the  separate  effects 
of  diffusion  and  convection  on  the  I.  N.  I. ;  in  contrast  to  the  far  wake,  convective 
effects  dominate  in  the  near  wake. 

To  augment  calculations  of  incipient  non -equilibrium  effects,  completely 
frozen  wake  flow  was  assumed  for  equivalent  cases,  and  the  validity  of  this 
assumption  was  then  tested  (equations  26  through  29).  The  results  for  the 
far  wake  are  shown  in  Figures  18  through  24.  Frozen  temperature  variation 
for  the  six  cases  along  the  centerline  is  given  in  Figure  18.  Figures  19  and 
20  show  the  variation  of  the  ratio  along' the  centerline  for  the 

assumptions  of  equilibrium  and  frozen  temperature,  respectively.  If  this 
rate  ratio  is  less  than  1  frozen  flow  is  justifiable.  These  results  show  that, 
considering  the  maximum  possible  temperature,  frozen  flow  does  appear 
justifiable  for  case  1  only.  Figure  21  shows  the  frozen  temperature  profile 
383^CjjA.  downstream  for  case  3;  Figure  22  shows  the  rates  Ej^  and  Ej^  and 
their  ratios  for  the  equilibrium  and  frozen  temperature  assumptions,  respec- 


tively.  These  results,  as  the  previous  one,  indicate  the  presence  of  a 
large  "intermediate"  flow  region  in  the  wake.  Figures  23  and  24  indicate 
the  effect  of  the  various  species  on  the  rates  for  the  same  profile,  the 
effects  of  0  and  N  appear  to  dominate  in  various  regions  of  the  wake. 

Finally,  Figures  25  through  27  show  similar  calculations  for  the 
four  near  wake  cases.  Frozen  temperatures,  and  the  rate  ratios  for  the 
equilibrium  and  frozen  temperature  assumptions  respectively  are  shown 
in  these  three  figures. 

SUMMARY 

Incipient  non-equilibrium  coefficients  were  found  to  be  a  practical 
and  versatile  research  tool  in  estimating  the  effects  of  finite  rate  chemical 
kinetics  in  complicated  hypersonic  wake  flow  models.  It  was  possible  to 
account  for  both  convective  and  diffusive  effects  in  a  multicomponent 
reacting  mixture.  The  incipient  non-equilibrium  coefficients  and  the 
application  of  chemical  kinetics  to  chemically-frozen  wake  calculations 
can  be  used  to  bracket  the  need  for  calculations  accounting  for  coupled 
fluid  mechanic  and  chemical  kinetic  effects  in  wakes.  The  results  indicate 
that  the  concept  of  "sudden  freezing"  may  be  useful  in  flow  regions  near 
the  body,  where  large  gradients  occur;  but  further  downstream  a  large 
intermediate  region  appears  where  neither  chemical  equilibrium  or 
chemical  "freezing"  can  be  justified. 
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NOMENCLATURE 


A 


C.,C. 
1  J 

AC. 


D_ 

Dt 

^chem 

E. 

1 

Z\E 

h 

H 

I.  N.I. 

n 

m 

M 

P 


Pr 


R 

Re 


reference  area  for  far  wake  drag 
mass  fraction  of  species 

th 

composition  deviation  (eqviation  8)  for  i  species 
drag  coefficient 

convective  derivative 

n+m 

chemical  energy  in  mixture  E  ,  =  )  E.  C. 

chem  L  11 

i=l 

chemical  energy  for  i^  species 
energy  deviation  (equation  11  and  12) 

energy  rates  due  to  reaction,  diffusion  respectively  (equations  27,  28) 

enthalpy 

total  enthalpy  H 

incipient  non- equilibrium  index  E/h  (equation  12) 

number  of  independent  species  in  mixture 

number  of  independent  conservation  laws  relating  species 

equation  (23) 

pressure 

equation  (9) 

Prandtl  number 

species  diffusional  flux  velocity 
far  wake  radius 
Reynolds  number 
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Sc 


Schmidt  number 


T 


u,  V 


w. 

1 


x.y 


6 

A 

e 


p 


chemical  and  fluid  mechanic  reference  times  respectively 

temperature 

equation  (9) 

axial,  radial  velocities  in  wake 

velocity  defect  in  far  wake  (equation  15) 

species  mass  generation  rates  due  to  chemical  reactions 

axial,  radial  coordinates  in  wake 

radius  of  near  wake 

equation  (23) 

expansion  parameter  (equation  2) 

equation  (23) 

viscosity 

density 

equation  (23) 


Subscripts 

0  condition  at  base  corner  in  inviscid  flow 

1,  2  first,  second  terms  in  6  expansion 
°o  free  stream 

b  base,  or  base  surface 

e  external  (inviscid)  flow 

eq  in  chemical  equilibrium 

Fr  frozen 

.  .  .th  .th 

1,  J  1  »  J  species 

ref  high  enthalpy  reference  condition  (equation  26) 
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INCIPIENT  NON- EQUILIBRIUM  INDEX 
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TEMPERATURE  PROFILES  IN  NEAR  WAKE  ENTHALPY  AND  CHEMICAL  ENERGY  IN  THE  NEAR  WAKE 
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VALIDITY  OF  CHEMICALLY  FROZEN  FLOW 
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The  effect  of  finite  chemical  reaction  rates  on  the 
chemical  structure  of  hypersonic  wakes  is  examined. 
The  technique  of  incipient  non-equilibrium  coefficients, 
based  on  linearized  chemical  kinetics  for  a  seven 
component  (Oj,  Nj,  O,  N,  NO,  NO\  ei  ten  reaction  model 
of  high  temperature  air,  is  utilized  in  the  analysis. 


The  far  wake  is  treated  as  an  asymptotic  region; 
the  near  wake  theory  considers  the  viscous  mixing  of  the 
boundary  layer  flow  with  the  inviscid  and  re-circulating 
flows  in  the  base  region,  and  the  subsequent  development 
of  this  viscous  core  into  the  wake  region  downst;  u.  m. 

The  results  show  where  non -equilibrium  effects  are 
important  in  the  wake  and  why  they  arise. 


The  incipient  non-equilibrium  results  are  then 
related  to  calculations  assuming  chemically  frozen  flow. 
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